Constraining Parity Violation in Gravity 
with Measurements of Neutron-Star Moments of Inertia 
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Neutron stars are sensitive laboratories for testing general relativity, especially when considering 
deviations where velocities are relativistic and gravitational fields are strong. One such deviation is 
described by dynamical, Chern-Simons modified gravity, where the Einstein-Hilbert action is modi- 
fied through the addition of the gravitational parity-violating Pontryagin density coupled to a field. 
This four-dimensional effective theory arises naturally both in perturbative and non-perturbative 
string theory, loop quantum gravity, and generic effective field theory expansions. We calculate 
here Chern-Simons modifications to the properties and gravitational fields of slowly spinning neu- 
tron stars. We find that the Chern-Simons correction affects only the gravitomagnetic sector of 
the metric to leading order, thus introducing modifications to the moment of inertia but not to the 
mass-radius relation. We show that an observational determination of the moment of inertia to an 
accuracy of 10%, as is expected from near-future observations of the double pulsar, will place a 
constraint on the Chern-Simons coupling constant of ^^^^ < 5 km, which is at least three-orders of 
magnitude stronger than the previous strongest bound. 

PACS numbers: 04.40.Dg,04.50.Kd,04.60.Rt,04.80.Cc 



I. INTRODUCTION 

Even though it has been almost a century since its 
original proposal, general relativity (GR) remains only 
marginally tested in the strong, dynamical regime, where 
velocities are relativistic and gravitational fields are 
strong. Great effort has gone into testing this theory 
in the solar system and with the binary pulsar (via the 
so-called parameterizedpost-Newtonian [l-Q and post- 
Keplerian frameworks ^^^^ ™^ cosmologi- 

cal settings (e.g., via the so-called parameterized post- 
Friedmannian framework With the imminent 

discovery of gravitational waves, new frameworks have 
been proposed (e.g., the parameterized post-Einstenian 
one |13| or the multipole expansion of black-hole space- 
times |14l - [l7| ) that will allow us to search for GR devi- 
ations with gravitational waves in the strong, dynamical 
regime. Until such observations become available, how- 
ever, one must rely on the next best tool to test GR in the 
strong field: observations of neutron stars and black holes 
in the electromagnetic spectrum (see e.g., Ref. [Tsl fioj). 

GR tests with neutron stars appear a priori to be not 
as clean as solar system or binary pulsar tests. When 
dealing with weak- field systems, we can usually employ 
the effacement principle [20, [2l[ to treat objects as effec- 
tive point particles moving under the influence of some 
force, which mathematically resembles Maxwell's electro- 
magnetism (the so-called gravitomagnetic analogy; see 
eg. [23]). In modified theories of gravity, however, the ef- 
facement principle and BirkhofF's theorem need not hold, 
but corrections usually arise at first order in weak-field 
perturbation theory (see e.g., for a review of internal 
structure effects in alternative theories). 



When dealing with neutron stars, we must inescapably 
account for the matter content of the star properly, as 
this can greatly infiuence the astrophysical observables 
to leading order. Our lack of detailed knowledge of their 
matter content is encoded in the equation of state (EOS), 
which relates the matter density to its pressure. Several 
EOSs have been studied in the literature (see eg. [l^HJ]), 
differing mostly in the nature of the strong interaction at 
very high densities and the presence of a softening com- 
ponent (such as hyperons, Bose condensates or quark 
matter), none of which we possess laboratory data for. 
Different possible EOSs thus lead to slightly different 
electromagnetic observables, which generically mimic al- 
ternative theory effects. Neutron star tests of GR must 
thus address possible degeneracies between the EOS and 
true GR deviations. 

Such a degeneracy between matter content and grav- 
itational theories is not new to GR tests. For example, 
the observed acceleration of the universe can either be at- 
tributed to a dark energy fluid or to modifications to the 
Einstein equations. This degeneracy in the explanation 
of observables is a consequence of the Strong Equivalence 
Principle, which essentially states that it is difficult to 
differentiate between a gravitational field and a matter 
distribution, since they are tied together by Einstein's 
equations. 

The only way to break this degeneracy is to possess 
at least two independent sets of neutron star observa- 
tions: one to determine the EOS as a function of the 
theory under consideration and another one to pin down 
the specific theory. Such two-observation tests have been 
dramatically successful when testing GR with binary pul- 
sars, where the degeneracies are associated with the un- 
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known masses of the pulsars [6|. New and exciting ob- 
servations are beginning to constrain the neutron star 
EOS [IB, [IB], suggesting that such two-observation tests 
will become possible in the near future. In this paper, we 
consider a particular EOS for neutron-star matter that is 
consistent with current observations. Our results, how- 
ever, depend weakly on our ignorance of the EOS. We 
show that the corrections introduced by the specific al- 
ternative theory considered couples strongly only to large 
density gradients, which occur only close to the neutron 
star surface, where our knowledge of the EOS is fairly 
robust. 

The construction of modified gravity theories is a non- 
trivial endevour that can be classified into two groups: 
infrared modifications and ultraviolet modifications (see 
eg-, [13, [13 for a review). The first group deals with cor- 
rections to the action that modify the weak-field predic- 
tions of GR. Examples of these include f{R) gravity |2S- 
l30l |. DGP gravi tv j31i] . Horava-Lifshitz gravity [32, l33| . 
and TeVeS [13, HH] . Usually, such modifications are in- 
troduced to propose explanations for the late-time ex- 
pansion of the universe or the anomalous rotation curves 
of galaxies. The second group corrects the action by in- 
troducing higher-order curvature terms, which, by con- 
struction, do not modify the leading-order predictions 
of GR in the weak-field. These theories do modify the 
strong-field regime of gravity, where neutron star obser- 
vations can place stringent constraints. Examples include 
Chern-Simons (CS) modified gravity [36*1, Gauss-Bonnet 
modified gravity scalar-tensor theories [38,- 4ft | . and 
certain /(i?) models [iU - lii ]. 

We concentrate here on tests of dynamical CS modi- 
fied gravity, which is currently only weakly constrained 
by binary pulsar observations 45| . This four-dimensional 
theory adds a certain contraction of two Riemann ten- 
sors and the Levi-Civita tensor to the Einstein-Hilbert 
action, coupled to a dynamical scalar field. It can be 
thought of as an effective theory (a low-energy approx- 
imation to some more fundamental theory) whose sig- 
nature is the modification of parity asymmetric gravita- 
tional solutions, such as the Kerr spacetime. It naturally 
arises as the low-energy approximation of many different 
fundamental theories. The latter could be string the- 
ory (i.e.. Type I, lib, Heterotic, etc.), where the inclu- 
sion of the CS term is inescapable in the perturbative 
secto r by the Green-Schwarz anomaly-cancelling mecha- 
nism [461 147| . In non-perturbative string theory, CS cor- 
rections also arise in the presence of D-instanton charges 
due to duality symmetries [48j . The more fundamen- 
tal theory could also be loop quantum gravity, where 
the CS correction has also been seen to generically arise 
when one of the fundamental constants of that theory is 
promoted to a field [U l49l - [5l| . Even without recurring 
to a specific fundamental theory, the CS term unavoid- 
ably arises in effective field theories as one of the only 
non-vanishing, second-order curvature corrections to the 
Einstein-Hilbert action in infiation [s^ . 

The non-dynamical version of CS-modified gravity 



(where the scalar field is non-evolving, but prescribed 
a priori) has been extensively studied and greatly con- 
strained in the past p8l [53l - l68| . Only recently, however, 
has there been a dedicated effort to analyze the dynam- 
ical theory [H, [gqI - ItH . In Ref. a small-coupHng 
approximation was developed to find the exterior gravi- 
tational field of a slowly rotating compact object. Per- 
haps rather surprisingly [72, 73], the solution was found 
to deviate from the Kerr spacetime because such a space- 
time sources a non-vanishing CS correction that couples 
to the spin-angular momentum of the black hole. Later, 
in [70|, the Strong Equivalence Principle was found to 
be satisfied in the dynamical theory and point particles 
were indeed found to follow geodesies of the background 
spacetime. 

In this paper, we calculate the CS-modified gravita- 
tional field inside neutron stars and relate this to possible 
observables that would allow us to constrain the theory. 
This analysis should be seen as complimentary to that 
of [i^, thus completing the gravitational field prescrip- 
tion everywhere in the spacetime of rotating bodies. Our 
analysis is also similar to that of [gB], except that here 
we focus on the dynamical theory instead of the non- 
dynamical one and we employ a more realistic, tabulated 
EOS. We shall see that if we employed the homogeneous 
stellar structure equation used in [63[ . the interior so- 
lution would not be CS-modified. Although in principle 
there will be CS modifications to the Dirac equation that 
describes the motion of fermions inside the star [gl] , we 
shall not consider those here. 

We find that the main CS correction affects the grav- 
itomagnetic sector of the metric to leading order, as was 
found for the exterior solution j45i] . Such a correction 
introduces modifications to the moment of inertia of the 
star, but not to the mass-radius relation to leading or- 
der. This implies that the usual degeneracy between al- 
ternative theory modifications and the EOS can be bro- 
ken: measurements of the mass-radius relation from sub- 
kilohertz neutron stars can be used to extract the EOS, 
while measurements of the moment of inertia can be used 
to test the theory, given the EOS. 

Such a degeneracy break is sensitive to the spin- 
frequency of the neutron star observed. If the neutron 
star is spinning faster than J7 ~ 1 kHz, then second-order 
corrections in the spin-frequency to the mass-radius re- 
lation are not necessarily negligible and the degeneracy 
is not broken. One might then wonder whether such 
rapidly rotating neutron stars exist in Nature. Up to 
date, most of the neutron stars observed spin slower than 
this limit (for a recent modeling of millisecond pulsars 
see [T^I)- However, over the past year there has been 
a breakthrough in millisecond pulsar searches through 
combined observations with the latest Fermi data. Such 
a breakthrough relies on millisecond pulsars generically 
being strong 7-ray emitters [tB]. If so, this pulsar popu- 
lation could be calibrated through the diffuse 7-ray back- 
ground they produce [77] (for an outline of such a new 
survey based on unidentified Fermi sources see fiBi , ) . The 
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tests we describe in this paper need not necessarily ap- 
ply to this population, as then the degeneracy between 
mass-radius relation and CS modifications could not be 
easily broken. 

Assuming that the degeneracy can be broken through 
observations of neutron stars with sub- kHz rotational fre- 
quencies, we find that a 10% accurate, moment of inertia 
measurement, as is expected in the near future with ob- 
servations of the double pulsar [tI, [t^, could place a 
constraint on the CS coupling constant of < 5 km. 
Such a constraint would be three-orders of magnitude 
stronger than the previous strongest bound [45j . 

Before proceeding with the remainder of the paper, we 
first comment on the relationship between neutron-star 
observations or tests and gravitational wave ones. In the 
past, these two programs could have been seen as com- 
peting approaches, when in reality they should be seen 
as complementary. Gravitational wave observations sam- 
ple the fully dynamical, propagating sector of a theory, 
while neutron star observations test the strong-field, but 
stationary sector. Moreover, simultaneous gravitational 
wave-neutron star observations, e.g. of a merging neu- 
tron star binary, could yield the most stringent bounds 
on alternative theories, as the combination of observa- 
tions would help break degeneracies. 

This paper is divided as follows: Sec. HIl defines the ba- 
sics of CS modified gravity; Sec. IIIII describes the metric 
parameterization and the stress-energy description of the 
neutron star; Sec. |IV] discusses CS modifications to neu- 
tron star models; Sec. IVl solves the CS modified equations 
of neutron star structure numerically and present some 
results; Sec. IVII computes some neutron star observables 
that are CS-modified and discusses possible future con- 
straints on the theory; Sec. I VIII concludes and points to 
future research. 

We use the following conventions in this paper: we 
work exclusively in four spacetime dimensions with sig- 
nature (—,-!-, d^; Latin letters (a, &, . . . , h) range 
over all spacetime indices; round and square brack- 
ets around indices denote symmetrization and anti- 
symmetrization respectively, namely T((ifc) = ^(T'ab-I-Tba) 
and T^ab] = 5 (Tab — Tfca); partial derivatives are some- 
times denoted by commas (e.g., 89/ dr = dr-O = 9^^), 
while radial derivatives are sometimes denoted with an 
overhead prime dritJ — lo' if w — ui{r)] the notation 
^(m,n) stands for a term of 0{m,n), which itself stands 
for terms of 0(e™) or 0{("); the Einstein summation 
convention is employed unless otherwise specified; finally, 
we use geometrized units where G = c = 1 and where 
Mq = 1.476 km stands for a solar mass. 



II. THE ABC OF CS MODIFIED GRAVITY 



A. Basic Equations 

The CS action is defined as 

S = SeH + ScS + 5*15 + Smat, 

where 
Scs 

*5'mat 



Kid x\/^gR, 
h 



- / d'^x^/^d *RR, 
4 Jv 



(1) 

(2) 
(3) 



1 



d^xV^ [g^' {Va^) (Vbz9) + 2V{^)] , (4) 

(5) 




Equation ([2]) is the standard Einstein-Hilbert term; 
Equation is the CS correction; Equation Q is the 
scalar-field term; Equation ([S|) describes additional mat- 
ter sources, with C^at the matter Lagrangian density. 
The quantity k^^ = IGttG is the coupling constant of 
GR, while a is the CS coupling constant. This formula- 
tion assumes that the CS scalar field 1!) is dimensionless, 
which then forces a to have units of length squared^. As 
usual, g is the determinant of the metric, Va is the co- 
variant derivative associated with the metric tensor gat, 
R is the Ricci scalar, *RR is the Pontryagin density 



KH — nn — H b n acd ; 



and the dual Riemann-tensor is 

1 



* rja cd . 



^cde f jja 
e R hef : 



(6) 



(7) 



with e'^^^^ the 4-dimensional Levi-Civita tensor with sign 
convention e0i23 _ -fl/^— where g is the determinant 
of the metric. 

The CS field equations are derived from the variation 
of the action with respect to the metric tensor and the 
CS couphng field: 



Rab 



— Cab — [Tab — TidabT 



dV a 



*R R. 



(8) 
(9) 



with Rab the Ricci tensor, □ = Va V the D'Alcmbertian 
operator, and Cab the C-tensor 

= (Vci9) eR^^ d + (VeVdi?) *R'ii'^by . (10) 

The total stress-energy tensor is 

(11) 



Ti T-imat I rp 

ab = J-ab + -'- ' 



1? 

ab-i 



In this section, we review some of the basics of CS 
modified gravity (see the recent review article ^ for 1 This formulation is equivalent to that of "i?] with the choice 
further details). /3 1. 
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where T"^^* stands for matter contributions, and T^^ is 
the scalar field contribution 

Tt, ^ (Va^) (V,^?)- iffab (Va^) {V^ ~ g^bV i^) . (12) 

We set the potential = at the scales of interest, 

with deviations possibly at the scale of supersymmetry 
breaking (see e.g., the arguments presented in [3q|). 

Taking the divergence of Eq. (jS]), the first term on the 
left-hand side vanishes by the Bianchi identities, while 
the second terms on both sides of this equation cancel 
each other because 

VaC''" = -i (V^z?) *RR = ^V'^T^,, (13) 

provided that Eq. ([9]) is satisfied. This then establishes 
that VaT^^^ = in dynamical CS modified gravity, thus 
forcing point particles to move on geodesies. 

Two versions of CS modified gravity exist: the non- 
dynamical version and the fully dynamical one. In the 
former, the scalar field is an a priori, prescribed function 
that lacks dynamics. With our choice of coupling nor- 
malization, this theory can be reproduced in the limit 
a — > oo, as then the scalar field evolution equation be- 
comes *RR = (the so-called Pontryagin constraint), 
which restricts the space of allowed solutions. We shall 
not consider this model here, as it has been heavily 
constrained by solar system and binary pulsar observa- 
tions [H, [13, m, [gI, [H, Hll, 113 • Instead, we shall con- 
centrate on the fully dynamical theory, whose equations 
of motion are given in Eqs. and 



B. Perturbative Constraints in 
Dynamical CS Modified Gravity 

The main idea of perturbative constraints (see e.g., [i^ . 
d^-d^l) is to perform a small-coupling expansion of the 
solution to a certain set of differential equations. From 
a mathematical standpoint, this is nothing but an appli- 
cation of asymptotic analysis or perturbation theory (see 
e.g., M)- 

In the context of CS modified gravity, perturba- 
tive constraints reduce to the small- coupling approxima- 
tion where one searches for solutions of Eqs. (|8]) 
and dHD in the limit C < 1, where ( is of Ola"^ / {C^ k)] and 
C is some characteristic scale that describes the physical 
system under consideration. For example, when search- 
ing for spinning black hole solutions j45j, £bh = AfsH 
with A/bh the black hole mass and C = / (kAI*). In 
the case of neutron stars, the natural length scale is the 
neutron star radius £ns — R and ( — /{kR"^). 

In addition to the use of the perturbative constraint 
scheme (i.e., the small-coupling approximation), we ad- 
ditionally focus on slowly-rotating solutions. We define 
slow- rotation by requiring that e <C 1, where e = 0{£ i?) 
and n is the angular velocity of the physical system. 



As before, since for neutron stars >Cns = R, then e = 
0{R n,), where is the star's angular velocity. 

The quantities e and ^ are book-keeping perturbative 
parameters associated with the slow-rotation and small- 
coupling approximations, respectively. We shall use them 
to remind ourselves of the order of the approximation. 
Note that ( ^ instead of ~ a because the evolution 
equation for z9 forces this field to be linear in a, and 
thus, the C-tensor becomes linear in a as well, leading to 
corrections to the field equations that scale as a^. 

Combining the small-coupling and slow-rotation ap- 
proximations leads to bivariate or two-parameter expan- 
sions. Schematically, all fields A are expanded as 

^ = E^"C"A™,")' (14) 

m.n 

where ^(m^„) is a term of C(e'", C"), which means a term 
of ©(e™) and ©(C")- For example, a term ^(i.i) is pro- 
portional to e C and thus of 0{e, C), while a term A(2,i) 
is proportional to ( and thus of 0{e'^,C,) [clearly A(2,i) 
is negligible relative to ^)]. The metric tensor can be 
expanded bivariately as 

yab — gab ^ ^SJab ^ '^^ab 

+ C^5r+^^9r+<^.^^+(^(3,3), (15) 

where 5^°'°'' is some background metric that satisfies the 

Einstein equations, while (<?ib'^\ 5'i),'°'') are first order 

perturbations and (<?ib'^\ si],'^'') are second-order 
perturbations. Given an exact solution to the Einstein 
equations, g^b^, its expansion in e must satisfy 

9^ = 9T+eg^'^ + e'9^"^+Oi^). (16) 

Thus, all terms g^^''^^ for all N are pieces of the slow- 
rotation expansion of a GR solution, while terms that 
are proportional to C represent CS corrections. 

The scalar field can also be expanded bivariately via 

^ab = ^T+C^T+e^T 

+ e^T+^'€b"^+^C^iY^+0{3,3). (17) 

Each piece in this expansion is determined by the evo- 
lution equation for the scalar field [Eq. ©], which in 
turn depends on the Pontryagin density. In [i^, [6^ . 
it was shown that for spherically symmetric spacetimes 
*RR = 0, thus forcing any physical CS scalar field to 
a constant and allowing us to set 'd'^^'^') = for all N. 
Angular momentum breaks the spherical symmetry and 
forces *R R ^ n ^ e to leading order, which then allows 
us to write 

^ab = e^i'f + eH'-^f + eC^iY^ + 0(3, 3). (18) 

We see that, to leading order, the C-tensor is linear in 
e, which then forces the leading order correction to the 
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metric to be proportional to (e (), i.e., a second-order 
correction in perturbation theory. The metric expansion 
simphfies to 



9ab 



GR 
9ab 



/■ (1,1) 



0(3,3). 



(19) 



Solving the modified field equations reduces simply to 
solving for 

The combination of the small-coupling and slow- 
rotation approximations establishes a well-defined iter- 
ation scheme. First, we solve the Einstein equations to 
linear order in e to obtain g^^'"'* -I- effi),'^'- Second, we 
use this metric to compute *i?i?(i^o) ^^nd from this solve 
Eq. dS]) for Third, we use this scalar field and the 

fifo'*^^ piece of the metric to compute the C-tensor, which 
when combined with the Einstein tensor G^),'^'' leads to 

differential equations for g^],'^''. Fourth, we solve these 
differential equations and iterate the scheme to higher 
order (see, e.g., Ref. [11] for more details on this boot- 
strapping technique). In this paper, we work to leading 
order only. 



III. METRIC PARAMETERIZATION AND THE 
STRESS-ENERGY TENSOR OF 
NEUTRON-STAR MATTER 

In this section, we consider the parameterization of the 
metric tensor and the description of the stress-energy ten- 
sor for neutron-star matter. We shall search for solutions 
to the field equations for a metric of the form [s^ 



d9^ + e^x sin^ 9 (dcj) - ujdtf 
(20) 

where (r, 9, (p) are (Boyer-Lindquist) spherical polar co- 
ordinates and (A,/3,7, x,w) are independent functions of 
radius and polar angle (r, 9) . The function uj is effec- 
tively the angular velocity of an observer falling in from 
infinity and, as such, is at least linearly proportional to 
the angular velocity of the rotating star. 

We decompose the metric functions in bivariate expan- 
sions as 



A 



A. 



(t)fi){r) A(2,o)(r,6') + eC A(i,i) (r, 61), 
P = /3(o,o)(^) + e'/3(2,o)(^^^) + eC/3(i,i)(^'^)' 
7 = 7(o,o)('^) + 7(2,0) (?-, 6*) + < 7(1,1) ('^'^)' 



X 



X(o,o) [r) + X(2,o) {r, 0) + eC X(i,i) ('', 



(21) 



To linear order in e, we can safely assume that all func- 
tions of 0(0, 0) and 0(1, 0) depend only on radius. No- 
tice that we have set (A(o,i), ^(o,i), 7(o,i), X(o,i)) = 
and (A(o,2),/3(o,2),7(o,2),X(o,2)) = as the C-tensor cor- 
rections start at 0(1,1), since there cannot be CS- 
deformations that are angular momentum independent 
due to parity invariance. Notice further that we have also 



set (a;(o,o),^^(2,o)) = and (A(i,o), /3(i,o), 7(i,o), X(i,o)) = 
because the GR gravitomagnetic and diagonal sectors are 
known to be series in odd and even powers of the an- 
gular velocity respectively in the slow-rotation approxi- 
mation [89]. Finally, we set 7(o,o) = 1 = X(o.o)i which 
amounts to a redefinition of the radial coordinate. 

In order to get a feel for how these functions behave, 
we first present the exterior gravitational field of a neu- 
tron star in GR. This exterior field is given by the slow- 
rotation expansion of the Kerr metric, where the small- 
rotation parameter e is now of 0{M a), with M the mass 
of the compact object and a its Kerr parameter. Notice 
that the Kerr parameter is related to the angular mo- 
mentum of the compact object via | J| — Ma, which then 
implies a = I fl/M, since by definition J = I fl, where 
/ is the moment of inertia (to be defined properly later) . 
To zeroth order in the spin, we find 



„2A\ 



'(0,0) 



1 



2M 



'(0,0) 



2M 



1> (e'")(o.o) = l' ^(0,0) =0. (22) 



'(0,0) ' \ Hofi) 
To linear order we find 



„2A\ 



'(1,0) 



- = (e^^ 



'(1,0) ' 



(1,0) 







=2x^ 



'(1,0) 



W(l,0) 



2Ma 



(23) 



To second order we find 



2Ma^ 



(2,0) 



(2,0) 



=27^ 



'(2,0) 



r.2 



cos 



(2,0) 



2M 
r 



2M 



cos^ 9 



2M 
r 



c^(2,o) = 0, (24) 



We model the stress-energy tensor of matter in the 
neutron star as that of a perfect fiuid 



(25) 



where p = p{r) and p = p{r) are the pressure and den- 
sity, while u"^ — 0, 0, r2M°) is the fluid's four-velocity 
that corresponds to a constant angular velocity Q, (we 
do not consider differentially rotating stars here). Nor- 
malization of this four- velocity requires Uau"^ = — 1, and 
hence 



+ ]^e^e-'^^'~'>'°''r^siTi^9{n-uj(i^Q)y 



(26) 



We now expand the stress-energy tensor order by order 
in e = RQ, ^ 1. To zeroth-order, the only non- vanishing 
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components are 



pe 



p e 



(27) 



(0,0)^mat _ 2 



{0,0)rpmat 
<P<t> 

sin^ e 



To first order in e, the non- vanishing terms are 

(l,0)ymat ^ _^2 gij^2 ^ (^f^ _ ^(1,0)^ ^ ^f^j ^ (38) 



while to second order in e we find 



(2,0)T^mat 



(2,0)Tnmat 
rr 

(2.0)Tnmat 



= 2A 



(2,0) P e 



2Ai' 



sin^ ( 



2 

(1:0) 



+ p^!'] +2pr2sin2 0(7(2,0) +X(2,o)) 



2pe2^ %o), 



g-2A(o,o) ^4 



(w(i'O) _f^^'(p^^)gi^4^^ (29) 



The cross terms of 0(1, 1) are 



-'-tt 

{l.l)rpmat 
rr 

-'-tip 
(l,l)^mat 



2 A(i,i) p e 



2Ai' 



/5(i,i), 
2 a 



r p sm 
2pr'^ sin2 6'(7(i,i) + X(i 



(30) 



The trace is always T™^* = g^-^T^^""^ = Sp-p up to ©(e^). 
Notice that the density and the pressure are assumed to 
be independent of f2 and a, and thus they carry no e or 
C dependence. 

The components of the stress-energy obtained in this 
section depend only on the perfect-fluid assumption of 
Eq. (|25p and the stationarity and axisymmetry assump- 
tion of the metric of Eq. pOI) . In general, for an arbitrary 
metric tensor, all components of the stress-energy will be 
non-vanishing, but here many components vanish due to 
the axisymmetry condition. Moreover, in the presence of 
electromagnetic fields, there will be additional contribu- 
tions that are not accounted for here. 

The EOS determines the relation between the den- 
sity and the pressure: p — p{p). We shall employ here 
the EPS EOS which predicts neutron-star masses 
smaller than 2Mq and radii (~ 10 km; see Fig. [T]) that 
are consistent with recent observations of X-ray bursters. 



IV. NEUTRON STARS IN DYNAMICAL CS 
MODIFIED GRAVITY 



In this section, we solve the modified field equations 
in the small-coupling and slow-rotation approximation 
(within the perturbative constraint framework). Each 
subsection considers the equations at each successive or- 
der in the perturbation theory ladder. 




8 10 12 
R [km] 

FIG. 1: The mass-radius relation for neutron stars calculated 
using the FPS EOS. Notice that for the FPS EOS, neutron 
stars are smaller than 2Mq and have radii ~ 10 km. 



A. Zeroth-Order Field Equations 

We first study the field equations to 0(0,0), i.e., the 
solutions that are independent of C and e. To this order, 
the C-tensor can be neglected, as it is at least linear in 
C and e. The field equations then reduce exactly to Ein- 
stein's equations (for a review of neutron star solution 
in GR see e.g. [H, [13, HH). The Einstein tensor pos- 
sesses only three independent components: (i,i), {r,r) 



and 



Using the definition m{r) = r{l — e p)/2 



the (i, t) and the (r, r) components reduce to the mass 
conservation equation and an equation for the A metric 
function, respectively 



2m', 



(0,0) 



a; 



(0,0) 



8np, 



Airpr^ 



r(r - 2m(o,o)) 



(31) 
(32) 



Hereafter, a prime stands for differentiation with respect 
to radius. We have explicitly included the index (0, 0) 
here to remind us that we are dealing with the 0(0,0) 
coefficients of the metric expansion. 

The system of differential equations becomes complete 
once we specify a third differential equation relating the 
rate of change of pressure to the density. We could use 
the (0, 9) component, but the equations simplify the most 
if we instead use the equations of stress-energy conserva- 
tion (which hold to all orders in dynamical CS modified 
gravity): \/^T™J^^ — 0. The r-component of this equation 
becomes 



^ fir 



^(0,0) {p + p)+p' = 0. 



(33) 



Combining this equation with Eq. (I32p we obtain the 
Tolman-Oppenheimer-Volkoff (TOY) equation 



^(0,0) 



Airpr^ 



r{r - 2m(o,o)) 



(34) 
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Given an EOS relating density and pressure, we can then 
solve the system of equations ([5T|)-([M1). 



B. Zeroth-Order Scalar-Field Evolution Equation 

Once we have obtained the metric functions A(o,o) 
/3(o,o)) we solve the evolution equation for the CS scalar 
field to zeroth order in e and ^. To this order, how- 
ever, the spacetime is completely spherically symmetric. 
In [6J|, it was shown that for such spacetimes, the Pon- 
tryagin density identically vanishes, i.e., *RR ^ 0, which 
then forces the evolution equation to become 

= . (35) 

Here again the (0, 0) index reminds us that we are work- 
ing to zeroth-order in both perturbation parameters. 

Since we are searching for a stationary and axisym- 
metric solution, can depend only on (r, 6) . With this 
assumption, the above differential equation becomes 



(0 



0) - e ■ff'^^y+iyy'e (^^(0,0) - P(o,o) 

(36) 



J_ „(0,0) cote (0,0) 



In fact, the above differential equation also holds to 
0(1,0) with the substitution ^ 

This differential equation only possesses ill-behaved so- 
lutions, i.e., solutions that lead to a scalar field with in- 
finite energy. Such an observation has in fact already 
been made when investigating the exterior gravitational 
field of slowly-rotating black holes . The only solution 
that is consistent with a scalar field with finite energy is 
^(0,0) _ const., which is the choice we make here. Of 
course, such a scalar field leads to no modification to the 
field equations as the C-tensor depends on i9 derivatives. 



C. First-Order Field Equations 

We now search for solutions of the modified field equa- 
tions that are either linear in e or linear in ^, i.e., solu- 
tions to 0(1, 0) or 0(0, 1). As already argued, however, 
there are no CS corrections independent of e, and thus 
'*^(o,i) — 0- As for corrections of 0(1,0), the C-tensor 
does not contribute, while the Einstein tensor contains 
only one equation, the (t, </>) component, which reduces 
in normal form to 



'^(1,0) ~^ ^^{1,0} + ^^{1,0) 



P + P 



I -2m 



(0,0) 



(37) 



Here, the dissipative and potential term respectively are 
4 



D 



V 



\'lr 



\o,o) ^ P[ofi) ) ' (38) 
2 
r 

2Afo,o) - 2A'('o,o) - 16^e2^(°^")p. (39) 



2/3(0,0) ^^'(0,0) ^ ~ ('^(0,0) ^(0,0)) 



We can simplify the above expressions using the zeroth 
order ones in Eqs. (PT|) - ([M)) . Doing so we find 



D 
V 



4 4:TTr^{p + p) 
r 



r - 2m(o,o) 



r — 2m 



(0,0) 



(P + P)- 



(40) 
(41) 



We consider some analytic solutions to this equations in 
the Appendix. Once more, once A, /3, p and p are known 
as functions of radius (from the previous rung of the per- 
turbative ladder), the above equations can be solved to 
find W(i 0) as a function of r. 



D. First-Order Scalar Field Evolution Equation 

We now consider the evolution equation for the scalar 
field to 0(1, 0) or 0(0, 1), i.e., we search for an evolution 
equation for this scalar field that is linear in e or linear in 
The left-hand side of the evolution equation is given 
by Eq. (jH]), as we showed in Sec. IIVBI The right-hand 
side of the evolution equation depends on the Pontryagin 
density, which to linear order is: 



*R R(i 



(1,0) 



r 

+ r 
+ A 



J-C0S6' e^^^w-o'-^C'") [l - e^''(°'«) 
(/3(o,o) ~ -^'(0,0)) + ('^'(0,0) ^ '^(o,o)/5(o,o) 



(0,0) 



"(i.o) 



(42) 



We see immediately that the Pontryagin density is of 
0(1, 0) and there is no 0(0, 1) piece, which immediately 
implies we can neglect t?(o.i). 

The complicated source for the evolution equation of 
the scalar field is greatly simplified once we make use of 
the lower-order equations. Inserting Eqs. (|3ip - (p4)) . the 
Pontryagin density simplifies to 



48 



3e-^<°.°) cos0(m(o,o) 



An 



-r p 



2m 



(0,0) 



(1,0)- 



(43) 



where w^-^ is to be understood as a known function 
obtained after solving Eq. (|37l) . For completeness, we 
present the full evolution equation for the scalar field 
[i.e., Eq. ([SJ to first order in e] below: 



2^(0,0,^(1,0) + ^(1.0),-2^<o,o, (^V^^_^^ _ p'^^^^^ + I 

-^(0,0) , 



(1,0) _^ C0ty^(i,o) 



12a 



HO,o) 



in , 



1 - 



2m 



(0,0) 



'(1,0) 



(44) 
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r [km] 



FIG. 2: The CS sourcing function m(o.o) — 47rr'^p/3 as a func- 
tion of radius in km. Close to the core, the function vanishes, 
while it reaches its maximum near the surface. The three 
curves correspond to stars with different central densities and 
hence with different masses. 



constant-density stars is the following. In '64'], it was 
shown that *RR = EatB"-'', where Eab and Bab are the 
electric and magnetic tidal tensors, which can be con- 
structed from the Weyl tensor. Such a relation allows 
for an analogy between CS modified gravity and elec- 
tromagnetic theory. In the latter, we know that the 
equations of motion for electromagnetic radiation satisfy 
wave equations in vacuum, but in the presence of matter, 
these equations acquire sources, e.g. DE cx Vp -I- J and 
OB cx V X J, where p and J are the charge and current 
densities respectively. Thus, in the presence of matter, 
the inhomogeneous solutions to the electric and magnetic 
fields need not be perpendicular to each other, leading to 
a non- vanishing E ■ B, and, by analogy, a no n- vanishing 
Pontryagin density. It is no surprise then that *RR is 
proportional to derivatives of the matter density and the 
current, the latter of which in the gravitational sector is 
proportional to derivatives of the cj metric perturbation. 



Before proceeding with the solutions, it is worth study- 
ing the structure of the Pontryagin constraint. The quan- 
tity in between parenthesis in Eq. (j43p is simply 



47r 



47r 



1(0,0) 



-r p 



[p - p) 



(45) 



where we have defined the mean density p = 
3m(o,o)/(47rr^). Clearly, for constant density stars, the 
above correction identically vanishes and the leading or- 
der term of the Pontryagin density is of This 
is indeed the case for the EOS employed in [6j], which 
would incorrectly suggest that for such an EOS there is 
no observable CS correction to neutron stars. Note here 
that there is no contribution to *R R because a term 
of this order is parity invariant, while only odd-powers of 
e can break this symmetry. 

For non-constant density EOSs, one expects the above 
density-dependent quantity to be close to zero near the 
neutron star core and reach its maximum value close to 
its surface, a behavior that is verified in Fig. [21 This 
figure shows the relevant term as a function of radius 
in units of kilometers. The three curves correspond to 
three different central densities that lead to different to- 
tal masses. Clearly, in the core (near zero radius), the 
Pontryagin density is nearly vanishing as the density is 
nearly constant, while it reaches its maximum near the 
surface of the star, which is close in magnitude to the 
total mass of the star. This behavior implies that the 
scalar field will be driven to non-vanishing values near 
the surface, which then in turn will lead to modifications 
to the gravitational field near the surface of the star. Due 
to this density dependence, the CS correction is not very 
sensitive to uncertainties in the nuclear physics, as the 
EOS near the neutron star surface is fairly robust. 

A physical reason behind the vanishing of *RR for 



E. Second-Order Field Equations 

The field equations to second order can be split into 
two sets: that of 0(1, 1) and that of 0(2,0) (as argued 
before, there is no correction to the metric of 0(0,2)). 
The latter reduces simply to the second-order GR equa- 
tions of structure, which we do not present here. This is 
a result of the C-tensor vanishing identically to 0(2, 0). 
The 0(1,1) equations contain CS modifications, as the 
C-tensor does not vanish in this case. 

The only surviving term of the C-tensor is 



(- 



Cu 



(1,1) 



-47rsin6' e^<°'°) f 1 - 



(46) 



+ ^ ) ^8^^(1,0) -ip-p) dre^ilfi) 



Since all quantities in this equation have already been 
computed, the C-tensor effectively acts as a new source 
for the field equations (a new effective matter term) . The 
accompanying component of the Einstein tensor to this 
order is given by 



G 



(1,1) _ 
tip 



— sin^ 0{1- 



2m 



(0,0) 



sm 



2r 1 - 



r 

2to(o 



(0,0) 



r 



- 2nr^ (p + p) 



Snpr'^ sin^ 6 i^(i.i) + ^^'^ ^ '^'^^ ^ de^{i,i) 



(47) 



The stress-energy tensor to this order has already been 
given in a previous section. The full (t, (f) component of 
the modified field equations to 0(1, 1) is then 
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2m 



(0,0) 



2r 1 



2m 



(0,0) 



Stt (p + p) 1^(1,1) = 47r 



Q,2 e^(o.") 



K smf 



r 
2m 



27rr^ (p + p) 



+ - (3 cot 6* deuji^i^i) 



(0,0) 



1/2 



^~V^ + ^ ) 5ei?(i,o) - (p - p) 9,.0Z?(i,o) 



(48) 



Given solutions for (A(o.o)j 't^(o,o)i '^(i. 0)1^1 P)? we can now 
solve the above equation for the CS correction to the 
metric 1^(1,1) • 

We can analyze separately the other components of the 
field equations. First, we note that all other components 
of the C-tensor identically vanish to this order. Moreover, 
the only non-vanishing components of the stress-energy 
tensor are the {r,r) and {<p,<j)) ones, which are 

linearly proportional to A(i_i), and (7(1,1), X(i,i)) 

respectively. The Einstein tensor also possesses non- 
vanishing components that lead to differential equations 
for these metric functions. However, since the C-tensor 
vanishes and there are no metric-independent matter 
source terms at this order, we can set these metric per- 
turbations to zero, i.e., (A(i4), ^(1,1), X(i,i), 7(1,1)) = 0. 



F. Boundary Conditions 
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0.7 - 

0.6 - 

"1 0.5- 

?^0.4- 

=^0 3- 
o 

0.2 - 
0.1 - 

0, 



0.2 0.4 0.6 0.{ 



M/M, 



1.2 1.4 1.6 



FIG. 3: The moment of inertia of neutron stars described by 
the FPS EOS, as a function of their masses. 



The differential equations presented in the previous 
subsections must be complemented by boundary con- 
ditions. The exterior solution for the metric outside a 
slowly-rotating star in dynamical CS modified gravity 
was found in 43. The line element was found to be 



dsi 



2 5 a a 

dsK + 7 4 

4 K 



5 a cos 9 
—a- 



— \- — ] sm"^ 9dtd(j), 



M 



2M 
r 



r 10 r2 



5r^ 



(49) 



where ds^ is the slow-rotation limit of the Kerr line el- 
ement. The solutions to the differential equations de- 
scribed in Sec. IIVEI must be guaranteed to approach 
the above solution in the limit r — > i?, where R is the 
neutron-star radius. 

In terms of the metric functions, the CS correction only 
affects the W(i,i) term, which becomes, for the exterior 
metric: 



5 a 



■'(1.1) 



12 M 
YV 



27 AP 
10~ 



(50) 



This relation is derived from the fact that 



,(1:0) 



— sin^ 9 a;(i q)- The metric components corresponding 



to the Kerr metric have already been given in Eqs. ((22)) - 
P4p . The boundary conditions are then simply applied 
by requiring that all metric functions approach their ex- 
terior modified Kerr counterparts at r = i?. 



V. NUMERICAL SOLUTIONS 

In this section we present results for the numerical so- 
lution to the modified field equations. We begin with the 
standard GR functions and proceed to solutions for the 
CS scalar field and the CS correction to the metric. 



A. General Relativistic Solutions 

We first solve numerically the GR equations of neutron 
star structure for the FPS EOS. We use standard numer- 
ical techniques to solve these equations: we reduce the 
system of equations ([3T|) . (|32l) . (|34|) and (|37| to first order 
and then employ a fourth-order, Runge-Kutta scheme to 
obtain a numerical solution. The boundary conditions for 
the GR solutions are implemented following the scheme 
of [i^, Is^l, which exploits the scale- free nature of the 
differential equations. 

We have already shown the mass-radius relation in 
Fig. [H From the gt^ metric component, we now com- 
pute the moment of inertia via 



Stt 




u^{r,9) 



dr . (51) 



Here gu and grr are the (i, t) and (r, r) components of 
the metric. We plot this quantity in Fig. [3] as a function 
of neutron-star mass. 
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B. First-Order CS Scalar Field Solution 



In order to solve Eq. (l44l) . we will first convert this 
partial differential equation to an ordinary one via sepa- 
ration of variables, i.e., 



,9(1^0) (r, ^) =^^„(r)P„ (cos ( 



(52) 



where P„ are the Legendre polynomials [62l and we have 
omitted the order symbols in for convenience. 
The 6'-dependent term in Eq. (|44)) is nothing but 



+ cote^^'/^ = (sine ^^l''^) . (53) 

' sm a \ ' / 

When we evaluate this with respect to Legendre polyno- 
mials, we find 
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FIG. 4: The absolute value of the dipolar component of the 
scalar field i? as a function of radius, for three stars with 
different masses and radii. 



I 



sm( 



■de (sin 6* P„,e) = -n{n + l)P„(cos6') 



(54) 



using the master Legendre polynomial equation. With 
this at hand, Eq. pij) becomes 



(0,0) 



COS 

(0,0) 1^ 



r 
An 



Pn 



(0,0) 



1 - 



2 m 



(0,0) 



(1,0) 



(55) 



Clearly, all 9 dependence disappears once we choose 
n = 1, such that Pi = cos 6* and 1?^^'°) = dcosO. This 
ordinary differential equation becomes 



12a X ( 47r o 



1 - 



2m 



(0,0) 



■^(1,0)' 

(56) 



which we can simplify even more using the lower-order 
TOV equations. The result is 



2m 



(0,0) 



-1 r 



— ( I 



2m 



(0,0) 



12a 



r V r 

^(°-«) ( m(o,o) - yrV 



2m(o,o)V'^' 
1 UJ 



(1,0)- 



(57) 



Equation ([5S)) is an ordinary differential equation (of 
dipole, n = 1, type) that must be solved for di. The 
boundary conditions for this function are the following: 
at the neutron star surface, cos 9 should equal Eq. (HHl) : 
at the center of the core, d'^ = 0, as there can be no CS 



correction when the spin angular momentum vanishes. 
This, of course, presents a numerical challenge, as the 
above differential equation is not scale free. We choose 
to use a shooting method, where we iterate over a variety 
of boundary conditions at the core for i9i until the surface 
boundary condition is satisfied. 

Figure [3] shows the absolute value of the CS scalar 
field in units of km~^ as a function of radius. In this 
plot, we have divided t9 by aQ to give it units of inverse 
kilometers, since a has units of squared length and Q of 
inverse length. The interior solution (r < R) is obtained 
by numerically solving Eq. (|57p . The exterior solution 
(r > R) is simply that of Eq. (gH). Note that the CS 
scalar overall increases in magnitude until the surface of 
the star, at which points it decays as the radius cubed. 
This occurs because as the radius increases toward the 
surface, the spin angular momentum increases until it 
saturates at r = i?. It is precisely this quantity (the spin 
angular momentum) that the CS scalar field couples to. 

The behavior of d is slightly different for the most 
massive star with M = I.SMq relative to the less mas- 
sive ones. The difference in behavior is rooted in the 
source term of the evolution equation for the scalar field 
[the right-hand side of Eq. (|57)) ]. This source is pro- 
portional to ui'f^i Qy which by definition is proportional 
to dlQYi/dr. Figure [3] shows that dlQj^/dr changes sign 
(ic. the slope of the curve changes sign) at around 
M = 1.4M0: dlcR/dr > for M < IAMq, while 
dlQu/dr < for M > 1.4M0. In turn, this implies that 
the source term in Eq. (|57p also changes sign. The change 
in sign of the source term forces d to start negative near 
the core, but quickly this quantity must switch sign so 
that the boundary condition in Eq. (|49)) is satisfied at the 
neutron star surface. This then explains why 9 behaves 
differently inside the neutron star for stellar configura- 
tions with masses larger or smaller than 1.4M0. 
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C. Solution to the Second-Order Field Equations chosen so that the left-hand side of Eq. p8|) simphfies to 



As before, we solve Eq. p8)) numerically. However, 
it is again convenient to convert this partial differential 
equation into an ordinary one via separation of variables, 
i.e., 



1 



— 9eP„(cos( 



(58) 



Such a functional dependence for the angular sector is 



3 cot 9 df, \ ^—rr I + a, 
sin( 



(n + 2) (n-1) 



sm 

( dePn 
V sin 61 



(59) 



after repeated application of Eq. (|54p . Using this relation 
and inserting Eq. ((5^ . Eq. (1151) becomes 



2to 



(0,0) 



sm( 



2r 1 



1m 



(0,0) 



Stt (p + p) a;„ 



sm ( 



-47r — e-^fo-f) 



2to 



(0,0) 



27rr^ ip + p) 

1/2 r 



aeP„ 



smf 



(n + 2) (n - 1) 



sm( 



P- P 



+ ^]^i-ip-p)^'i 



(60) 



where, in the right-hand side, only the n — 1 term 
survives in the expansion. Once more, for the left- 
hand side of this equation to be independent, we must 
choose n = 1 in the w*^^'^^ expansion, which then ren- 
ders The above differential equation then 
simplifies to 



( ^ 2m(o,o) 



2r 1 



2m 



(0,0) 



27rr^ (p -I- p)] - Stt (p + p) = 47r — e^*" 



o \ 1/2 



(61) 
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FIG. 5: The absolute value of the CS correction to the metric 
coefficient lo\ as a function of radius inside and outside a 
neutron star, for stars with three different masses and radii. 



Given that ?9i is known, the above equation is simply an 
ordinary differential equation for wi. 

The boundary conditions for Eq. (pT|) are the follow- 
ing: at the neutron star surface, Lb\ is given by Eq. (I50p : 
at the center of the core, = 0, as, once more, there 
can be no CS correction when the spin vanishes. As be- 
fore, we employ a shooting method, where we iterate over 
boundary conditions at the core for uj\ until the surface 
boundary condition is satisfied. 

In Fig. [5] we plot the absolute value of the CS correc- 
tion to the gravitomagnetic metric component in units of 
km~^ as a function of radius for three stars with different 
masses and radii. The interior solution presents a rather 
constant behavior for the two less massive stars, while 
the exterior solution always decays with the sixth power 
of the radius, as expected. 



VI. CONNECTION TO OBSERVABLES 

In this section, we connect the above CS modifications 
to observable quantities, in order to investigate whether 
neutron star observations can constrain the CS coupling 
parameter. For simplicity, let us define ^ = / k, such 
that C, — ^/ R^. A constraint on the CS coupling a is thus 
a constraint on ^. It is customary to quote constraints 
on cx a^/^, as this quantity has units of length. 

The mass continuity equation and the equation of hy- 
drostatic equilibrium are CS modified to subleading or- 
der, i.e., 0(2, 1), because the uj expansion (see Eq. (|2T|) ) 
cannot have a 1^(0.1) piece and any correction to the 
equation of hydrostatic equilibrium is proportional to 
~ ^lifl) + 2a;(i^o)'^(i,i) + 2^^(1,0)^^(2,0) • This implies 
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that CS modifications to the mass-radius relation scale 
{^/R'^)iM/R){mf, which is clearly 
suppressed by a factor of ^ 1, where the factor of 

must be present for dimensional consistency^. If the 
mass radius relation where measured to an accuracy of 
<5m-r ^ 10%, we could then place a constraint of order 



,1/4 



< 55 km 



<^M-R 



R 



0.1 



10 km 

1/4 



0.025 



1/2 



\RQ/c 

R 1.6Mq 
10 km M 



1/4 



(62) 



< 



which derives simply from requiring that (5m-r 
w(i,i)a;(i,o)/(Af/i?) ~ i^/R^) {Rnf . We have scaled the 
previous expression by i? = 10 km, M — I.GMq and 
Q/{2tt) — 600 Hz, corresponding to the typical radius 
and mass and the largest spin frequency of neutron stars 
bursters for which we have measurements psj . Of course, 
for more slowly rotating pulsars, for example if ~ 45 Hz 
(ri ~ 0.36 Hz), corresponding to the A (B) component 
of the double binary pulsar (79| , then the bound worsens 

dramatically: CntR < 220 km i^Htn < 2400 km). 

The weakness of the CS correction to the mass-radius 
relation of slowly-rotating stars is in fact an advantage. 
For sufficiently slowly-rotating stars, for example, for 
stars whose spin frequency is $7 ^ 10'^ Hz, corrections 
of 0{e^) to the mass-radius relation can be ignored. In 
such cases, the degeneracy between EOS and CS correc- 
tions is broken. That is, given a measurement of the 
mass-radius relation for a sub-kHz neutron star, the CS 
correction can be ignored, and the mass-radius relation 
can be used to infer the EOS. 

A stronger constraint arises from the measurement of 
the moment of inertia of a neutron star, as is suggested 
by the modification to the gravitomagnetic sector of the 
metric [see Eq. (1511) ]. Since the metric function uj is 
CS modified, we infer that / is also modified. We then 
decompose / = /qr + Ics j where the latter is given by 
Eq. (ISTj) with the substitution (1 — oj/H.) — ^ 
and {gtt,grr) given by their GR values. Since the mo- 
ment of inertia scales by definition as l/il, the CS modi- 
fication is first-order in e. We can then infer that roughly 
I ^ Igr (l + S./R'^), and assuming that / has been mea- 
sured to an accuracy Sj (and found to agree with the GR 
result), we could infer the following constraint on ^: 



< 5 km 



R 



10 km 



0.1 



1/4 



(63) 



Here we have assumed again a fiducial neutron-star ra- 
dius of 10 km and the nominal accuracy for the error in 
the measurement of the moment of inertia of 10% (78l . [79l |. 



^ This relation can also be obtained by substituting back into 
'^{i,0)'^(i,l) the surface values required by the boundary con- 
ditions. 
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FIG. 6: The ratio between the CS correction to the moment 
of inertia and the GR value as a function of the CS coupling 
parameter ^^'^^ in units of kilometers. The thick horizontal 
line corresponds to a measurement of the moment of inertia 
with an accuracy of 10%. 



As expected, this constraint is independent on the rota- 
tion rate of the neutron star. 

Another way to estimate the magnitude of the con- 
straint is to consider directly the lo metric perturbation. 
Since the moment of inertia depends linearly on oj, a mea- 
surement of / to an accuracy Si is effectively a measure- 
ment of UJ to the same accuracy. Given then a measure- 
ment w, we can infer oj {1 + S) = a;(i.o) (l + i^(i.i)/i^(i.o)) j 
which leads to the constraint S < g), or simply 



< 5 km 



R 



3/4 



M 



1/4 



1l 

0.1 



1/4 



10 km/ V1.4Af0 

(64) 

where we have here used Eqs. ([25]) and to model 
the mass and radius dependence. As expected, this con- 
straint is consistent with the one shown in Eq. (|63|) . 

In order to obtain more specific constraints on the CS 
coupling, we solved for the CS correction to the moment 
of inertia numerically. In Fig. [B]we plot the ratio of the 
CS correction to the moment of inertia to the GR value as 
a function of ^^Z'' for three stars with different masses and 
radii (the lowest-mass case is shown for illustrative rea- 
sons, as no astrophysical process is known that leads to 
such light neutron stars). We remind the reader that we 
have previously defined ^ = c? j k. A 10% measurement 
implies a constraint of roughly t^l^ < (5.5, 6.5, 7) km for 
stars with final mass M = (0.8, 1.4, 1.8)Af0. These con- 
straints are of the same order of magnitude as the rough 
estimates presented above. Also notice that the expan- 
sion parameter ( — ^/R* is of order 0(10^^) if we satu- 
rate these constraints, which justifies the small-coupling 
approximation. We should emphasize that, unlike con- 
straints derivable from measurements of the mass-radius 
relation, these constraints are independent of the spin 
frequency of the neutron star. 

The possible bounds that one could achieve with mo- 
ment of inertia measurements is to be contrasted from 
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that obtained from observations of the orbital evolution 
of the double binary pulsar CbinLry ~ km [i^. The 
moment-of-inertia bound will be much stronger because, 
in the binary pulsar case, the CS coupling parameter cou- 
ples to the binary separation (roughly 10^ km), which is 
much larger than the neutron star radius. 



VII. CONCLUSIONS 

We have studied the structure of neutron stars in dy- 
namical CS gravity. We have found that to leading or- 
der, only the gravitomagnetic sector of the metric is CS 
corrected. This implies that the mass-radius relation is 
only modified to second order, while the moment of in- 
ertia is modified to leading order. We have derived the 
CS modified equations of stellar structure and we have 
solved them numerically. These solutions complete the 
prescription of the gravitational field inside and outside 
of stars. We estimate that a measurement of the moment 
of inertia to 10% can lead to a constraint on the theory 
that is at least three orders of magnitude stronger than 
previous constraints. 

As is typical with tests of alternative theories with 
neutron star observations, the observables are degener- 
ate with the EOS, which leads to a degeneracy between 
alternative theory modifications and our knowledge of 
nuclear physics. In the case of CS gravity, however, the 
mass-radius relation can be used to extract the EOS, as 
neither the mass nor the radius are CS corrected to lead- 
ing order. This is simply a consequence of the fact that 
CS gravity couples to leading-order to the spin angular 
momentum and not to the mass density. 

Future observations that lead to a measurement of the 
moment of inertia of a neutron star would make the test 
presented here viable. The resulting constraint holds the 
promise of being one of the strongest achievable with any 
type of astrophysical observation, at least in the static 
sector of the theory. This analysis could then be com- 
plemented by gravitational wave tests, which sample CS 
corrections in the radiative sector. 
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Appendix A: Some Analytic Solutions to the 
Gravitomagnetic Sector of the Metric 



Equation ([37l) governs the behavior of the gravitomag- 
netic metric perturbation W(i_o). Unfortunately, there is 
no known closed-form, analytic solution for any realistic 
EOS. However, for some particularly simple EOSs, such 
a solution can in fact be found. 

One example is to set p = — p, corresponding to a 
gravastar configuration j94. l95j. i.e., a gravitational vac- 
uum star, whose interior solution is diffeomorphic to De 
Sitter spacetime. In this case, the dissipative function 
becomes D = 4/r, the potential V = and the right- 
hand side of Eq. (1571) vanishes. We can then solve for the 
gravitomagnetic potential to find 



^^(1,0) = ci + 



C2 



(Al) 



where ci.2 are two constants of integration. Finally, we 
can also solve the equation of hydrostatic equilibrium to 
obtain p = C3, where C3 is a constant of integration. 

Perhaps a more realistic configuration is one with van- 
ishin g pr essure and constant density p ~ pc- In this case, 
Ref. [631 has found that the GR solution is 



6 / r \ 2 
5 VR) 



(A2) 



assuming the fiuid is rigidly rotating with angular veloc- 
ity vl. 

A third alternative is to allow p + p = {r — 



2m 



(0,0) 



)/(7rr3). In this case, = and F = -U 



while the right-hand side of Eq. ((37)) becomes equal to 
— 16il/r^. The gravitomagnetic component becomes 



(A3) 



where pi^2 — 1/2 ± V65 and ci,2 are constant of integra- 
tion. 
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